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Relation between Events in the Millimeter-wave Core and Gamma-ray Outbursts 
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Analysis of comprehensive monitoring of 34 7-ray bright quasars, BL Lac objects, and radio 
galaxies reveals a close connection between events in the millimeter-wave emission imaged with the 
VLBA at 43 GHz and flares at 7-ray and lower frequencies. Roughly 2/3 of the flares are coincident 
with the appearance of a new superluminal knot and/or a flare in the millimeter- wave "core" located 
parsecs from the central engine. This presents a theoretical challenge to explain how the 7-ray flux 
can often be variable on intra-day time-scales. Possible answers to this include very narrow opening 
angles of the jet, small volume filling factors of the highest energy electrons, chaotic magnetic fields, 
and turbulent velocity fields relative to the mean jet flow. 



INTRODUCTION 



A pre-requisite for understanding the violently vari- 
able nonthermal emission in blazars is to identify the 
location of the high-energy emission region. Many 
arguments implicate a relativistic jet pointing almost 
along the line of sight, but there has been a heated de- 
bate among theorists as to whether most 7-ray flares 
take place within the sub-parsec broad emission-line 
region (BLR) [e.g.,|3l| or within and downstream of 
the "core" seen in millimeter- wave very long baseline 
interferometry (VLBI) images [e.g., |23|. There has 
been less controversy among observers, who have con- 
sistently found too many coincidences of 7-ray out- 
bursts with radio events on parsec scales to support 
the notion that 7-ray flares occur on sub-parsec scales 
[3 113) EH- Nevertheless, the statistics regarding the 
fraction of correlated events relative to the number 
expected by chance, as well as the relative timing of 
the maxima of 7-ray and radio flares, have stimulated 
skepticism among many of those who need BLR pho- 
tons to generate 7-rays [e.g.,|3l| and steepen the 7-ray 
spectra above 2 GeV from opacity to pair produc- 
tion [29.]. 

The all-sky survey mode of the Fermi Large Area 
Telescope (LAT), plus the concerted effort by many re- 
search groups to provide commensurate data at other 
wavebands, are combining to provide new data sets 
that are sufficiently rich to decide between these two 
scenarios. Here we report the results of over three 
years of roughly monthly monitoring of 34 7-ray bright 
active galactic nuclei with the Very Long Baseline Ar- 
ray (VLBA) at 43 GHz (wavelength of 7 mm). At this 
frequency, the VLBA, an array consisting of ten an- 
tennas scattered around the continental US as well as 
the islands of St. Croix and Hawaii, provides an an- 
gular resolution as fine as 0.1 milliarcseconds (mas). 
This translates to 0.6 pc at a redshift of 0.5. Since 
the jet on these scales is usually not highly opaque at 
43 GHz, the VLBA gives us a clear view into regions 
quite close to the central engine. 



II. OBSERVATIONS 

A given VLBA session extends over 24 hours, dur- 
ing which time we observe 30 objects spread out in 
right ascension. (Some of the sources with relatively 
slow proper motions are observed once every second 
month). The integration time on each source ranges 
from about 30 to 60 minutes, split into segments of 4- 
8 minutes so that a wide range of spatial frequencies 
(baseline vectors of pairs of antennas in wavelength 
units) is sampled. After calibration in AIPS, we edit 
the data in Difmap to remove bad points when an 
antenna experienced an obvious problem. This is fol- 
lowed by iterative creation of an image via an inverse 
Fourier transform of the visibility data, followed by 
application of the CLEAN algorithm to construct a 
pattern of flux density within a grid of pixels — the 
model — and then self-calibration with heavy con- 
straints to match the data better with the model. Af- 
ter several iterations, we form an image by convolving 
the model with an elliptical Gaussian function that 
approximates the central portion of the point-spread 
function ("beam"). After generating a high-quality 
preliminary image, we transfer the CLEAN model 
back into AIPS for self-calibration of the phases, with 
the routine CALIB adjusting the left and right cir- 
cularly polarized data independently. This removes 
residual differential phase errors. Since a rather good 
image of the source serves as the basis for the ad- 
justment, the binning time of the phases is allowed 
to be short (6 seconds). We then transfer the cor- 
rected data back into Difmap for final imaging and 
self-calibration. Further steps are required to cali- 
brate the polarization; these are described elsewhere 

[a. 

The images resulting from our VLBA monitor- 
ing program are available at the website of the 
Boston University blazar research group at URL 
www.bu.edu/blazars/VLBAproject.html, along with 
the calibrated uv data and CLEAN component mod- 
els. The images and data are available to other re- 
searchers who are interested in individual sources or in 
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applying their own analysis and interpretation. Publi- 
cations using the results should acknowledge that our 
data have been incorporated into the study and cite 
the above URL. 



III. LOCATION OF GAMMA-RAY FLARES 
IN THE JET 

Figure [T] displays a sample of our VLBA data set, 
a five-epoch sequence of VLBA images of the quasar 
3C 273. The range of epochs is selected to coincide 
with the first stage of a major, multi-flare 7-ray out- 
burst. The bright, compact feature at the northeast 
end, the "core," is used as a reference point since ab- 
solute positional information is absent owing to un- 
known delays of the wave front by the Earth's atmo- 
sphere. The core has the characteristics of a stand- 
ing conical shock 0, HO] ■ Figure [2] shows multi- 
waveband light curves over a three-year period, with 
the dates marked when a new superluminal knot was 
coincident with the core. Four knots were ejected 
during the multi-flare 7-ray outburst starting around 
RJD 5050. Just as significantly, no new bright knots 
appeared during the 450 days following the end of the 
outburst, when the 7-ray emission stayed at a low 
level. 

We have examined our VLBA images of all 34 
sources for both variability of the emission features 
and the appearance of superluminal knots ejected 
from the time when Fermi began science operations in 
August 2008 through August 2011. Our prehminary 
analysis finds that 43 7-ray flares were simultaneous 
within errors with the appearance of a new superlu- 
minal knot or a major outburst in the core at 7 mm. 
This compares with only 13 cases (in 11 blazars) in 
which either a 7-ray flare occurred without a corre- 
sponding mm- wave event, or a mm- wave flare or su- 
perluminal knot ejection had no 7-ray counterpart. 
(Four of these blazars had other 7-ray flares that did 
correspond to a mm- wave event.) In addition, flve 
blazars that were quiescent over the three years at 
7-ray energies were also quiescent at 7 mm. Radio 
associations with eight 7-ray flares were less obvious, 
with several months separating isolated radio and 7- 
ray events that only occurred 1-2 times in three years. 
In at least two of these, AO 0235-^164 and OJ287, the 
7-ray flares coincided with superluminal knots cross- 
ing quasi-stationary (i.e., slowing shifting) features in 
the jet. Such coincidences, combined with a robust 
statistical analysis, confirms the association of the 7- 
ray and millimeter- wave events 0, [4*| . 

The evidence is therefore clear that most 7-ray 
flares occur near or in the mm-wave core. But not 
all of them do so. An example where the evidence in- 
dicates that some of the flares took place upstream of 
the core is PKS 1510-089 0,[22|. An apparently con- 
tinuous rotation of the optical polarization vector by 
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FIG. 1: VLBA images at 43 GHz of the parsec-scale jet of 
3C 273 during the early stages of the series of 7-ray flares 
that started around June 20, 2009 (RJD=5003). Two new 
superluminal knots, K4 and K5, appeared over 2 months. 
The linear polarization (false color) in the jet reached max- 
immum on RJD=5091, which coincided with the highest 
peak in the 7-ray light curve. Total intensity contours are 
in factors of 2 starting at 0.4% of the peak of 3.68Jy/beam. 
The maximum in polarized intensity is 0.115 Jy/beam. 
The restoring beam is an elliptical Gaussian with FWHM 
dimensions of 0.38 x 0.15 mas, with maximum elongation 
along position angle —10°. 
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FIG. 2: Gamma-ray and X-ray light curves of the quasar 
3C 273 from August 2008, when regular 7-ray monitoring 
by Fermi began, to mid-2011. The vertical arrows in the 
bottom panel indicate times when the brightness centroids 
of new superluminal radio knots passed the brightness cen- 
troid of the core on 7 mm VLBA images (see Fig. 1). Note 
how the rate of appearance of knots became high during 
the main 7-ray activity. 



720° over 50 days in 2009 implies that a single emis- 
sion feature was responsible for the highly variable 
emission during this time interval. During these 50 
days, there were several optical and 7-ray flares with 
different ratios of the fluxes of the two wavebands. 
An extremely sharp, high-amplitude (brightest level 
observed since 1948) optical flare and equally sharp 7- 
ray flare (of unknown maximum flux, since the blazar 
was not observed much of that day owing to an au- 
tomatic re-point of the Fermi LAT detector) erupted 
while a new superluminal knot was passing through 
the core. These events are easiest to explain either if 
there are local sources of seed photons along the jet, 
or if the seed photon field is highly time- variable [I^l • 
The discovery that most of the highly variable 7-ray 
activity arises in — or sometimes even downstream 
[1) 01 — of the 43 GHz core leads to a very profound 
realization: the VLBA at millimeter wavelengths is 
able to image the j-ray emitting region! It is there- 
fore an extraordinarily powerful probe of the physi- 
cal conditions that cause 7-ray flares. Although our 
enthusiasm is tempered a bit by the weak, if any, re- 
lationship between TeV flares and events in the jet 
in high-spectral-peak BL Lac objects [23], we are en- 
couraged by the connection between such TeV flares 
and new superluminal knots in 4C21.35 and BL Lac 
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FIG. 3: Fermi LAT 7-ray (blue) and 43 GHz VLBI "core" 
(red) light curves of ten of the quasars in the observing 
program. Yellow rectangles: times (± uncertainties) when 
a new superluminal knot passed the core. Thin blue arrows 
indicate 7-ray upper limits. 



IV. STRONGLY VARIABLE 
HIGH-FREQUENCY RADIATION ON 
PARSEC SCALES 



There has been considerable resistance to the con- 
clusion that the majority of 7-ray flares occur parsecs 
downstream from the central engine [e.g.,[3l|. This is 
understandable, since the 7-ray and optical, and often 
X-ray, flux is usually highly variable [e.g., 111. How can 
intra-day time-scales of flux changes [e.g.,l3| be consis- 
tent with a parsec-scale source? Even with high levels 
of relativistic beaming — the Doppler factors 6 of the 
most extreme blazars can exceed 50 [13| — the time- 
scale of variability of a spherical source of radius Rpc 
parsecs at redshift z is limited by light-travel argu- 
ments to longer than tvar,min ^ 20i?pc(l + z){S/50)~^ 
days. 

The answer to this question must be that the vol- 
ume involved in the highest-energy emission is much 
smaller than the corresponding length scale of the jet. 
The first argument in favor of this is that the most 
highly relativistic jets — those of the brightest 7-ray 
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FIG. 4: As in Figure [S] Fermi LAT 7-ray (blue) and 43 
GHz VLBI "core" (red) light curves of another ten of the 
quasars in the observing program. 



blazars — are extremely narrow. For example, careful 
analysis of 43 GHz VLBA data allowed Jorstad et al. 
[13} to derive S ^ 25 and an opening semi- angle of 0.8° 
for 3C 454.3 at z = 0.86. If the 43 GHz core lies - 10 



pc from the vertex of the jet cone then at this 
point the jet has a cross-sectional radius of only ~ 0.14 
pc. This translates to a lower limit of ivar.min 10 
days. So, a time-scale of ^ 0.5 days for flux decline 
from the peak of the November 2010 mega-flare in 3C 
454.3 3 requires that only a fraction ~ [0.06((5/25)]^ 
of the cross-section is involved in the 7-ray emission. 
If this extraordinary event was stimulated by a par- 
ticularly high Doppler factor, 5 ^ 50, then a small 
region covering roughly 1/60 of the jet cross-section 
was responsible for the 7-ray emission. 

Following a preliminary model that showed promise 
[13], the lead author has developed a turbulent ex- 
treme multi-zone (TEMZ) model [over 3000 zones; 
that attempts to explain how such events can occur. 
The overall variability is caused by fluctuations in the 
magnetic field strength and particle density at the jet 
base, perhaps from instabilities in the accretion disk's 
magnetic field, described via a red-noise power spec- 
tral density. The jet plasma is turbulent, with many 
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FIG. 5: Fermi LAT 7-ray (blue) and 43 GHz VLBI "core" 
(red) light curves of ten of the BL Lac objects in the ob- 
serving program. Yellow rectangles: times (± uncertain- 
ties) when a new superluminal knot passed the core. Thin 
blue arrows indicate 7-ray upper limits. 



cells — each with uniform magnetic field — across the 
jet cross-section. Electrons are accelerated to high en- 
ergies, with a power-law distribution, as the turbulent 
plasma crosses a standing conical shock. The number 
of turbulent cells that radiate synchrotron emission 
at a given frequency can be estimated from the level 
of linear polarization. For example, a mean optical 
polarization of ^ 10% can result from ~ 60 cells if 
the magnetic field orientation is random from cell to 
cell. Since electrons with similar energies should pro- 
duce both optical synchrotron and several hundred 
MeV inverse Compton emission, this should also cor- 
respond to the number of cells involved in the 7-ray 
production. If one of these cells has a particularly 
high density of electrons, then it can cause a sudden 
7-ray/optical flare as it crosses the conical shock if the 
acceleration of electrons is particularly efficient. This 
might require a particular orientation of the magnetic 
field relative to the shock front [s^. The variability 
can be accentuated if a turbulent velocity field is su- 
perposed on the mean flow of plasma down the jet 
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V. DISCUSSION 

The TEMZ model conforms with a number of in- 
dications that the old picture of strong^^^ transversely 
oriented shocks propagating down jets 9, can de- 
scribe only some flux outbursts and superluminally 
moving emission features. The position angle of lin- 
ear polarization is often oblique or nearly perpendic- 
ular to the jet axis [e.g., knots Kl and K3 in Figure 
[1] see also |25|, rather than parallel to the jet axis 
as expected for such shocks. In addition, increases 
in flux are often simultaneous with passage of mov- 
ing knots through stationary bright features that are 
consistent with standing shocks [3, [l^ . While propa- 
gating oblique shocks can explain the polarization po- 
sition angles [13], the red- noise power spectra of the 
flux variations at radio, optical. X-ray, and 7-ray fre- 
quencies [l|, seem more consistent with stochastic 
processes than with a small number of singular events 
such as strong shocks. There is also ample evidence 
that turbulence plays a major role in the variations: 
the optical polarization of blazars is typically of order 
10% and it is highly variable, characteristics expected 
of synchrotron radiation in a turbulent magnetic field 

The data lead us to propose that most of the vari- 
ability of flux and polarization variability of blazars 
and appearance of superluminal knots is caused by 
random fluctuations in the magnetic field (magnitude 
and direction) and electron density and energy dis- 
tribution stemming from noise processes. (Some of 
these fluctuations in energy density might steepen to 
form shocks as they propagate down the jet.) Flares 
across a large portion of the electromagnetic spectrum 
occur as turbulent cells in a region of higher than nor- 
mal energy density cross standing shocks. Extremely 
rapid flares at the highest energies (optical to X-ray 
for synchrotron radiation and 7-ray for inverse Comp- 
ton scattering) can arise when a particular turbulent 
cell is highly efficient at accelerating electrons to very 
high energies, perhaps because its magnetic held is 
favorably oriented relative to the shock front. 

VI. CONCLUSIONS 

Comparisons of 7-ray with radio light curves and 
the appearance of new superluminal knots have long 



suggested that in blazars the two opposite ends of 
the electromagnetic spectrum are intimately related. 
Now, combined monitoring with Fermi and the VLB A 
conflrms this close association. This conclusion poses 
a great challenge to theorists to flgure out how vari- 
ability on time-scales as short as hours (or minutes in 
some TeV flares) can occur parsecs from the central 
engine. The lead author is pursuing a model involv- 
ing turbulent plasma flowing across standing shocks 
that shows promise in reproducing these and other 
observed characteristics of blazar emission. The most 
extreme cases may require explosive events such as 
magnetic reconnection [8|. The outcome should be a 
more comprehensive theory of the physics of relativis- 
tic jets. 

So, the spirits of Fermi and Jansky appear to be 
working together to create a glorious celestial show. 
This is only fltting, since blazars dominate both the 
7-ray and millimeter sky at high Galactic latitudes. 
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